Far-from-equilibrium superconductor in fluctuational regime 
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We derive Ginzburg-Landau-like action for two-dimensional disordered superconductor under far- 
from-equilibrium conditions in a fluctuational regime. Then, utilizing it, we calculate fluctuation 
induced density of states, Maki- Thomson and Aslamazov-Larkin type contributions to the in-plane 
electrical conductivity. We apply our approach to thin superconducting film sandwiched between a 
gate and a substrate that have different temperatures and different electrochemical potentials. 
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I. INTRODUCTION 

Most of processes in physics and in technological realm 
occur under far from the equilibrium (FFE) conditions. 
At the same time, the theories of nonequilibrium behav- 
ior were mostly restricted to small deviation from the 
equilibrium. A marked progress in approaches to quan- 
titative description of FFE physics is related to Keldysh 
technique-based methodstl. However, mostly FFE sys- 
tems have been studied within the nonequilibtium form 
of the mean-field theory, sec, for examplc,afl. In the 
present paper, we address the important question of what 
happens to the conventional second-order phase transi- 
tion if the system under consideration is driven out of 
equilibrium. In equilibrium when close enough to the 
second order phase transition the mean field theory does 
not hokLand the physics starts to be governed by fluctu- 
ations.liEl In our work, we construct a theory of fluctua- 
tions under FFE conditions. As the exemplary system, 
we consider a superconductor FFE in the fluctuational 
regime. 

Our first notion is that while in clean three dimensional 
conventional superconductors the fluctuations are impor- 
tant only in a very narrow region around the supercon- 
ducting transition line (usually within the ~ 10 _12 K tem- 
perature range), in high temperature, low-dimensional 
and organic superconductors the fluctuation region is 
much wider. In particular, as early as in 1968, Aslama- 
zov and Larkin, and, independently, Maki showed that 
in disordered thin superconducting films the width of the 
fluctuation region, which is determined by the sheet re- 
sistance, grows noticeably as compared to that of bulk 
superconductors. En3 Moreover, it was demonstrated that 
not only thermodynamic but also dynamic characteristics 
of the low-dimensional systems are strongly influenced by 
the fluctuations close to an equilibrium, see e.g. Refs. 0,0 
for a review. 

A quantitative approach to nonequilibrium fluctuation 
superconductivity was recently formulated in Refs. foUO. 
Building on this approach, we develop further our origi- 
nal Keldysh technique enabling us to find nonequilibrium 
fluctuation contributions to the electrical conductivity of 



a superconductor above the (nonequilibrium) supercon- 
ducting transition. We show that by measuring the fluc- 
tuation corrections one can infer the parameters of the 
nonequilibrium state of the superconductor from the ex- 
perimental data. Further, while in the equilibrium the 
lifetime of the fluctuation induced Cooper pairs is deter- 
mined by therdifferencc T — T c , where T c is the critical 
tcmpcraturcElEl, we find that in FFE conditions it is con- 
trolled by the parameters of the nonequilibrium density 
matrix of the system. For example, for a thin supercon- 
ducting film sandwiched between the gate and the sub- 
strate, see Fig. EL these parameters are the temperatures 
of the gate and of the substrate, and the gate voltage V G . 
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FIG. 1: Thin superconducting film sandwiched between the 
substrate and the gate. The substrate temperature is Tf, the 
gate temperature is T2 and the gate voltage is Va- Changing 
Ti, T2 and Va one can tune the nonequilibrium quasiparticle 
distribution in the film. 



The paper is organized as follows. In Sec. [ffl we intro- 
duce the model. In Sec. [II we derive the nonequlibrium 



Ginzburg-Landau (GL) like action, while in Sec. [V 



calculate fluctuation induced corrections to the electrical 
conductivity of the normal metal. In Sec. M we focus on 
the specific experimental realization shown in Fig. |l| and 
analyze the previously derived results for this concrete 
setup. In Sec. VI, we summarize our results and discuss 
their further applications. Technical details are relegated 
to Appendix A. 



II. KELDYSH SIGMA-MODEL ACTION 

The Keldysh partition function of a superconductor 
with the Hamiltonian H in the coherent state basis is 
defined as: 



Z=M \ D^D^exp{iS[ij,ij}}, 



where 



s®M 



dt 



(ipaidtipa) - H[%l),ll)] 



(1) 



(2) 



and A/" is a normalization constant. Hereafter h = c = 
ks = 1. Here C is the Keldysh contour and a =t,4- is the 
spin variable. 

The Hamiltonian has the form H = Ho + -Hint, where 
the single-particle Hamiltonian 



H o = / ip a 



(V - ieA) 
2m 



U dh 



1>a. (3) 



Here A, <f> and C/dis are vector, scalar and disorder poten- 
tials and e is the electron charge; the tensor summation 
over the spin indices is implied. The interaction Hamil- 
tonian describes the electron-electron interaction in the 
Cooper channel, 



Him = -- J i>^i>\.i>i^, 



(4) 



where p(r) = ^ Q V'a( r )V ; a( r ) is the local electron density 
and the superconductive coupling constant A > 0. The 
disorder potential is assumed to be Gaussian distributed 
with the correlator 



1 



<E/dis(r)t/di s (r')) = <5(r - r') 



1~KVT 



(5) 



Averaging over disorder and carrying out the standard 
decoupling in the four-fcrmion terms in_jthe action via 
the Stratonovich-Hubbard fields {Q, A},EI and integrat- 
ing out the degrees of freedom with the energies higher 
than 1/t, where r _1 is the elastic scattering rate, we ar- 
rive at the Keldysh nonlinear cr-model partition function 



Z= I D[Q, A] exp{»S[Q,c, A*]}, (6) 

S[Q^A,c} = S A + S^ + S Q . (7) 

The nonlinear cr-model action 5* consists of three parts: 



Sa 



2A 



Tr[A]rY& K ], 



Scj, =—Th:[<j> K Y(j) K ], 



ITTV , 



Sq =^Tr[D{d r Q K ) 2 - 4:Zd t Q K 



(8) 
(9) 

(10) 



where Q 2 — 1 . Here D is the diffusion coefficient and v is 
the bare single particle density of states at the Fermi level 
per one spin projection. The action (mj holds while the ef- 
fective temperatures (see Sec. |TV| ) that follow from it are 
much smaller than 1/t. The check mark above the field 
variables indicates that they are defined in the space that 
is the tensor product of the Keldysh and Nambu spaces. 
The former and later are spanned by the Pauli matrices 
&i and fj, % € {0, x, y, z}, respectively So, Y = a x (g> fo, 
S = <7o (8> t z . Multiplication in time space is implicitly as- 
sumed, and "Tr" includes an integration over real space. 
The subscript JC denotes the gauge transformed fields: 



4>K.=<j)- d t JC, 

Ak = A + V/C, 
t = [k cl & + k q a x 



TO- 



A and <p are defined in same way as /C. Also 

A = [A c; <T + A q a x ] (g> t+ - h.c, 



(11) 

(12) 
(13) 

(14) 
(15) 



Qx. is defined in the same way. The quantum (q) and 
classical (cl) components are defined as half-sum and 
half-difference of the field values at the lower and upper 
branches of the Keldysh time-contour. The field A cl be- 
comes the superconducting order parameter on the mean- 
field (saddle-point) level, while the saddle point equation 
for Q produces the Usadcl quasiclassical equations, where 
Q plays the role of the quasiclassical Greens function. 
The covariant spatial derivative is given by 



9 t Qk = VpQk; - ie[EA/c, Q K }. 



III. GINZBURG-LANDAU ACTION 



(16) 



First, let us consider simplified situation ignoring the 
interactions. Then, the metallic saddle point of Eq.(ffl), 



A = UA U \ A = (7 z <8>f z , 

V U — Ot-t'+OTo 



Ft A*) 



F?Av) 
F t \,(r) 



(17) 
(18) 

(19) 



By setting quantum components of electromagnetic po- 
tentials to zero, the equation for F^ t , (r) reads: 

D(y 2 F e ' h T 2ie[A c l, VF e / h ] T ie[VA c l,F e ' h ] 



- e -[A c l [At F e ' h ]]) - tF e ' h + F e / h K 
T ie[<t>tF e ' h ]=Q. 



(20) 



After the Wigner transformation, we can map F e ' h 
to quasiparticle electron (hole) distribution functions: 



r e/h. 



e e/h 



Fj (r, t) = 1 — 2/ e ' (r, t). In the leading order with 
respect to time and energy partial derivatives of external 
potentials, the kinetic equations become 

D [ V 2 F e/h T 2e(d t A c ^)(d e VF^ h ) 

+ \e 2 {d t A%f{d 2 F*/ h )\ -d t F e J h 



T ed t {DVA% + ^)d e F^ h = . 



(21) 



Note that one has to take into account adequate bound- 
ary conditions. 

Having specified the metallic saddle point solution, we 
consider the masslcss. fluctuations around it. They can 
be parameterized aai3 



Mr) = 
W = 



-W(r)/2 



A(r) e 



W(r)/2 



WT + — W T- WqTq + W Z T Z 
WqTq + W Z T Z WT + — W*T_ 



W = UWU~\ (22) 
(23) 



such that WX + AW = 0. Here we introduced four real 
fields w" t , (r) , w" t , (r) with a = 0, z representing diffuson 
degrees of freedom and two complex fields w-u> ( r ) ; wtt' (r) 
for Cooperon degrees freedom. 

Since the main goal of this section is the derivation 
of the nonequilibrium extension of GL-like action in the 
normal state but very close to the transition, in what 
follows we will concentrate on Cooperon degrees of free- 
dom only. By plugging Q^, given by Eq. (|22j), in the 
action Eq. (|7p and expanding it up to the second order 
in Cooperons we get 



Sq — S u 

■KV 

dm? = i— Tr 



vA 



w| ]t ,(r)C-}(r)w t ,, t (r)' , 

(-A cl * - F h A q *)w + (A cl 



+ (A cl * - F e A q *)w + (-A cl + F h A q )w* 
Here w M /(r) = (w t ,t>{r),w ttt ,(r)) T ancH 

C- t )(v) = -dt + d v - le [tf£( T ,t) - ^(r,f)] 
-D[V-ieA^(v,t)-ieAi(v,t f )} 2 , 

C- t )(r) = d t -d t ,+ie [^(r,i) - <^(r,i')] 
-D[V-ieA^t)-ieA^{v,t')} 2 . 



(24) 
(25) 
F e A q )w* 
(26) 



(27) 

(28) 
(29) 



Note that there are additional terms in the action (24), 
that are not presented here for the sake of brevity. They 
contain quantum components of electromagnetic poten- 
tials. Since quantum components of fields are the aux- 
iliary source fields usually used to calculate observables 
by the appropriate differentiation of the Keldysh action, 
we omit the mentioned terms here. They are not impor- 
tant for the derivation of the GL action. We will discuss 



these terms in the next section, when calculating differ- 
ent corrections to the Drude conductivity that arise due 
to superconducting fluctuations. 

Now, the Cooperon degrees of freedom can be inte- 
grated out from the Keldysh partition function (pf). In 
this way the GL action is generated: 

iS G , [A ci ,A?] = -±(Sl A ) iSa2 + iS A - (30) 

In what follows, we consider stationary distribution func- 
tions ( i.e. Ft,t> = Ft-f) and stationary electromagnetic 
fields. Also, we assume that they slowly vary in space 
as compared to the fluctuating order parameter A(r, £). 
Then, we obtain (see App. [A]) 



5 GL [A C/ ,A"] =2^Tr AJ^r,*)!" 1 A K (r, t) 



(31) 



where A = (A cl , A q ) T . The structure of the fluctuation 
propagator L is characteristic for bosons: 



_ L A 

T- 1 T- 1 



(32) 



The indexes R, A and K denote the retarded, advanced 
and Keldysh propagators, respectively. The cl-cl com- 
ponent of L is zero. This is expected, since for A q = 
the field on the upper and lower branches of the Keldysh 
contour is the same, and therefore the corresponding ac- 
tions cancel each other, resulting in S GL (A d , 0) = 0. The 
Keldysh propagator is responsible for the Gaussian noise 
term in the time-dependent GL (TDGL) equation, and 
also, it carries information about distribution of electrons 
and holes. The details of the derivation, as well as general 
formulas for the fluctuating propagators (that are valid 
if the system is far from the transition) are presented in 
the App. El Close to the transition into superconducting 
state, the propagators become 



^oy> 



Lk = * 2 ' l 

T- 1 - — 

R / A ~ 8T e 



(33) 



TeF R - 

£)}■ 



T MT e F? + £>( V - 2zeA^) 2 T d t T 2ie<$ 



x l±i 



(34) 



Here F, 



:/h 



F 



:/h 



cl denotes the gauge invariant distri- 
c±ei f > /c 
bution function, while other parameters appearing in the 

retarded and advanced propagators are the functionals 

of-Ff= (F^-F± e )/2: 



T- 



AF 



R 



de 



e=0 
R 



9 r F H — F R 

n- 1 = --/de e „ ° . 



cp 



1 



(35) 
(36) 

(37) 



They are strongly drive dependent, as will be demon- 
strated in Sec. [V| The symbol f denotes the principal 
value of the integral. Note that the existence of effec- 
tive temperature T e , does not imply the local equilibrium 
form of excitation distribution functions. 

The nonequilibrium GL relaxation rate is defined as 



"■^cp^ e 



+u D Pf - tanh 
de 



2T C J 



4z cp §F «, 



(38) 



where ujd is the Debye energy. The GL relaxation 
rate represents the inverse lifetime of Cooper pairs, and, 
therefore, it vanishes at the transition to the supercon- 
ducting state. In an equilibrium, the density matrix 
is parametrized by the temperature, and the condition 
r^" L = tells us that the transition occurs at T = T c . 
In a general nonequilibrium case, additional parameters 
may appear, e. g. voltage drop and temperatures of the 
thermal baths that are in contact with the system. Then, 
the condition t~^ = defines the phase transition surface 
in the parameter space. We find that any distribution 
function F e at the phase transition surface satisfies: 



+un F e R 
de 



tanh I Tpp- 



■21 



*P" 







0. 



(39) 



Tuning the external electromagnetic fields and/or tem- 
perature of thermal baths, one can control the distance 
from the phase transition surface. We point out that 
the theory presented above is valid only for the systems 
(in the normal state) close to the transition, i. e. when 
(^GiTe) -1 < 1. The difference between the Eqs.(pi| 
and the work in Ref.JTol is the appearance of O 



38|) 
and 



presence of terms ~ Fq [results of Rcf.|l(]| are restricted 
to the nonequilibrium states generated by the voltage 
between the leads or gates in a certain way]. In general, 
fl^ 1 is nonzero, as discussed in Sec. JV|. 

In an equilibrium the following relations hold: F e ' h = 
tanh (^p\, T e = T, fi = 0, z cp = 1 and we reproduce the 
standard GL action. It is given by Eq. (pw , but with the 
propagators: 



(l; 



)K = l-T, 



(40) 



(L-^r/a = ^ [-Kir 1 + DV 2 T dt] , (41) 



(rSD =-(T-T c 



(42) 



for lo <C T. We find that in a general case, this relation 
is violated for a system out of the equilibrium. 

Next, we derive the TDGL equation. After the term 
~ |A 9 | 2 in Eq. (|3l| ) is decoupled by introduction of the 
Hubbard-Stratonovich field £(r, t), one differentiates (pL) 
with respect to A 9 * and obtains: 



^•X^Ag + C-O, 

7T 



(44) 



16, 



(C(r,i)C(r',i')) = — Tl ( 1 - F^) S(t r')S(t t'). 

(45) 



The field C, has a meaning of the Gaussian noise. The 
superconducting order parameter changes under a gauge 



transformation as A^ 



A c 'e 2iefcC '. However, the TDGL 



equation for the order parameter is gauge invariant, since 
the presence of the gauge field k cl in L^ 1 (through A^ 
and 4%) compensate the change of the order parameter 
phase. 

Using the Keldysh formalism, derivation of the GL 
close-to-equilibrium action was performed in RcfJ9|. 
There was utilized a special gauge where D'VA/c + <\>k = 
and found that the scalar potential does not appear in 
TDGL equation for the order parameter. In RcfJfl it was 
used F e = tanh(e/(2T)). By looking at kinetic equation 
( pl| ) it becomes obvious that mentioned gauge may in- 
deed simplify the calculation. However, in the presence 
of external potentials tanh(e/(2T)) is not the solution of 
the kinetic equation even in this gauge, due to boundary 
conditions. As a result, the scalar electromagnetic po- 
tential does appear in the retarded propagator, (34), and 
therefore also appears in the TDGL equation for the or- 
der parameter, (ff4|), contrary to the statements in Reffi 

Although the higher order terms in the superconduct- 
ing order parameter in the GL action are unimportant 
for our analysis of superconducting fluctuations, we will 
state them for completeness. We focus only on local 
fourth order terms of the form A 9 *A c/ |A d | 2 and find: 



S 



A 1 



i I <l\ttY(t) 
lo 



x Tr 



A£M)A£(r,i)|A£(M)r +c.c, (46) 



Y(r) 



K 

de F C R 



_^ 2ire + i0 



In the equilibrium the action reduces to 



-7K(3)TV A£ (r, t) Ag (r, t) \ A£ (r, t) \' + c.c. 



(47) 

J ^A^ = 

/(4tt 2 T 2 ). 



In an equilibrium, the Keldysh fluctuation propaga- 
tor satisfies the relation following from the fluctuation- 
dissipation theorem 



(L-^k = coth (£-) [(L-^nCq.w) - (L£) A (q,u>)] 



\2TJ 



(43) 



IV. CORRECTIONS TO ELECTRICAL 
CONDUCTIVITY 

In this section, we derive corrections to the Drude con- 
ductivity arising in thin films due to superconducting 
fluctuations. The system is close to the transition, but 
under FFE conditions. 



First, we briefly explain the origin of different fluctu- 
ation induced contributions. Fluctuation Cooper pairs 
carry charge and directly contribute to the electri- 
cal conductivity, determining Aslamazov-Larkin (AL) 
contributionQ'Ll. Since quasi-particles are involved in 
Cooper pairing, effectively the number of carriers partici- 
pating to the single-electron charge transfer is decreased. 
The fluctuation pseudo-gap opens at the Fermi level in 
the single-particle spectrum, and results in diminishing of 
the Drude conductivity for the so-called density of states 
(DOS) correction!! The third contribution to conduc- 
tivity, the (anomalous) Maki-Thompsonmla (MT) correc- 
tion is of purely quantum origin. It arises due to coherent 
scattering of the electrons, forming a fluctuation Coper 
pair, on the elastic impurities. As a result, the diffusion 
coefficient changes and therefore the Drude conductivity 
is influenced. 

We are interested in linear response to the in-plane 
electric field, while the system is strongly driven out of 
equilibrium due to contacts with thermal baths at dif- 
ferent temperatures between which it is sandwiched, or 
due to an electric field perpendicular to the plane. We 
focus on the regime where fluctuations can be treated 
pcrturbatively, meaning that below derived corrections 
are small compared to the Drude conductivity, a D . We 
take into account terms linear in the Ginzburg number 
Gi = (isDdf)^ 1 , i.e. 5<j/<j d ~ Gi, where df is film thick- 
ness. The in-plane dc conductivity is given by 



— lim 

2 n->o 



f 



S 2 Z 



n*Ag_(n)«*A« (-n) 



where Z is the Kcldysh partition function, Eq. (|fj|). The 
terms in the action, that are linear or quadratic in the 
vector potential, arise from the first term of Eq. (|I0|) and 



read: 



iS A = -^-Tr{[W,{VW)}ZA K } 



(49) 



iS^ = 



irvDe 2 



Tr{[HA K , A][5Ajc, A] + 2BA K hEA K KW 2 



2SA K KW 



SA K Ajy}, 



(50) 



respectively. The first term of Eq. (^) docs not contain 
superconducting fluctuations and therefore after differ- 
entiation over classical and quantum components of the 
vector potential and averaging over Q and A fluctuations 
(see Eq. ( 48 ) ) , it g ives the Drude conductivity. The other 
two terms in (p0[), produce the MT and DOS correction 
to the conductivity. The part of the action linear in the 
vector potential, Sa, gives the AL correction. The Drude 
conductivity is 



■jrvDe 2 
4 

-nvDe 2 
2 
2vDe 2 . 



lim < — Tr (A e a x An +e + & x K e A- 



n+e. 



lim 
0->-o 



Tr (Ff 



pe , ph 



F h \ 



o 



(51) 

In Eq. (ffq), all the quantum components of vector poten- 
tial arc set to zero after differentiation. Therefore, after 
differentiation with respect to A cl and A q , in calculation 
of a we can use the action ( |24|) derived in the previous 



, (48) section. We assume stationary situation (F i 



e/h 



1 t-V r 



The in-plane classical components of the vector potential 
can be set to zero, since we are interested in the liner re- 
sponse to the in-plane electric field. Then, we find the 
saddle point configuration of w for the action (|24|) and 
using the GL-like- action (|lj), we obtain: 



J 



((wi,2(q)u>3 i4 (-q)»Q,A = — S ei 

2i 

((wi, 2 (q)w* 4 (-q))) Q , A = — S ei 

1% 

«u?i, 2 (q)i4,4(-q)»Q,A = —S ei - 

1% 

((wi, 2 (q)w3. 4 (-q)»Q,A = —6, 



—L K L A ,1-2Lr : 1-2 + ^3 ifi,l-2 + FiLa.1-2 

e2 ' i4 ~ €3 {Dc i 2 -i{e 1 +e 2 )}{Dc i 2 -i{^ + e i )} : 

—L K La,1-2Lr,1-2 — F 2 La,1-2 — -PT^it',,1-2 

e2 ' i4 ~ €3 {Dq 2 + i{ ei + e 2 )} {Dq 2 + i(e 3 + e 4 )} : 
L K La,\-iLr,\-'i + F 2 La.i-2 — F 3 1 Lr^i^2 

" £3 {£>q2 + i( ei+C2 )}{£> q 2_ i ( e3 + e4 )} ' 



ei-C2,e4-e3 



L K La,\-iLr,\-i + F%Lr.i-2 — F{La.\~2 



(52) 
(53) 
(54) 
(55) 



r 



Here the angular brackets denote averaging over fluctu- 
ations of Q,A c '/«. F\ 



e/h 



r e/h 



and w.j 



i, j = 1 . . .4. Also, L A /R,i-j = [L R ) A (q, e* - e,) 



Now, we can proceed to calculation of different cor- 
rections to the Drude conductivity. Wc start with DOS 
correction: 



Sa DOS = — 



vDe< 



8- 



((^3-F3 e )((^3,4(qX )3 (-q) + ^3,4(q)< 3 (-q))) QiA 

+ (^3^)((w3,4(q)^4, 3 (-q) + W* 3A {(l)w 4 .3(-q)))Q.A 

- Fi l d n ((w 3 A(q)wl +n3+n (-q)))\ n=Q - Fida((w3 A (q)w i+ ci,3+a(-q))) 
+ F 3 e an((w3,4(q)w2_ ,3„ n (-q)))| n=0 +F 3 ' l an((w5 i4 (q)w4-n : 3-n(-q)}}| 0= 



n=oJ' 



(56) 



r 



where/ = Jdq/(27r) 2 J"_ de3de4. Note that terms 

that in the limit O — >• behave as 1/17 will be canceled 
out with similar terms from the other corrections (MT 
and AL), and that is why they are omitted here. The 
main contribution close to the transition reads 



6a z 



De 2 

47T 



111! 



LJ? (d e F?\L$\ 2 + d e F?\L*\ 2 ) 



{Dq 2 - i(2e - uj)Y 



(57) 



where L± = La(<1, ±w). We obtain the DOS correction 
in a thin film 



5a L 



7e 2 C(3) T e T cp 



n 4 df T 2 \ t 

" "/ I DOS \"gl 



T„ 



log ^r h ( 58 ) 



up to logarithmic accuracy. Here df is the film thickness 
and new characteristic temperatures are 



T, 



cp 



T e z cp (l - F^), 



T 2 

-'DOS 



7C(3) 



Re 



d e F L 

de- ^-—t; 

(e + iO) 2 



(59) 
(60) 



Fl = (F e e + F^)/2. The important contribution in 
Eq. ( p7|) comes from small momenta, and therefore we 
safely cut the momentum integration at the upper limit 
D( lliax ~ T e- The main contribution in the DOS, MT 
and AL corrections comes from frequencies uj — 2e<j)+eo < 
t~2 • Here and in the following we assume that the sys- 
tem is close to the transition, such that characteristic en- 
ergy scales of F e 
e = -4T e F R 



e/h 



are much greater than t~^ and |eo|/2; 
T~^T e /Vl. If the system is far from the 
transition and these conditions are not satisfied, then one 
can start calculation from Eq. (p9) and use general fluc- 
tuation propagators (see App. JA|) that are not restricted 
to low frequencies. 



Next, we focus on the MT correction to the conductiv- 
ity. It is given by: 
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The main contribution close to the transition is: 
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(62) 
(63) 



where we cut-off the infrared divergency in the momen- 
tum integration by introduction of .-the finite dephas- 
ing time Dq min ~ t$, t gl < t^.eI There are many 
phase-breaking sources, such as the electron scattering on 
phonons or paramagnetic impurities, or superconducting 
fluctuations.u The nonequilibrium conditions may affect 
also the equilibrium phase breaking time. However, we 
leave this problem for future studies. By treating energy 
e as a complex number in (|62j), we obtain the first term 
in Eq. ( J63] ) from the poles of the integrand determine d by 
zeros of the denominator. The second term in Eq. (|63] ) 
comes from the poles of distribution functions F e / h . Note 
that the first term is positive, while the second one is neg- 
ative. Then, the DOS correction is effectively doubled, 
although the first term in 5a mT is the dominant one close 
to the transition. 



Next we calculate the AL correction to the conductiv- 



ity: 
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(64) 



(65) 



We find that close to the transition the main contribution 
assumes the form 






d 
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where the new characteristic temperature is given by 
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t a \ ^ 2 Vr a 2 nr a T b J' 



= Im 
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(67) 
(68) 

(69) 



Performing the remaining integration over q, u, we ob- 
tain in the quasi-two-dimensional case: 



da 






(70) 



We conclude that all the fluctuation-induced corrections 
to the conductivity behave differently as a function of 
t gl . For a thin film close to the superconducting tran- 
sition, the MT is the most important one. Moreover, 
each correction is parameterized by different combina- 
tion of the effective temperatures: T e , T cp , T DOS and T AL . 
These temperatures are strongly drive-dependent, as will 
be shown in the next section when considering a con- 

? e/h 



In this section we have focused on the derivation of the 
fluctuation conductivity corrections close to the transi- 
tion. They are the most pronounced in that region, but 
nevertheless they can be still significant also far from the 
transition^ Then, the derived GL-like theory is not ap- 
plicable. One has to take into account high-frequencies 
and short-wave contributions in fluctuating propagators. 
However, this can be done within above developed ap- 
proach. Namely, the expressions for the corrections given 
by Eqs. (|5^,f3l],M) are valid also far from the transition. 
Then, in Eqs. ( p2f|55|) one has to use the general expres- 
sions for the fluctuating propagators, that are given in 
the App. |a|. 



V. EXPERIMENTAL REALIZATION 

In this section we propose an experimental setup where 
our predictions could be tested. The setup is shown in 
Fig. l|, where the superconducting film is sandwiched be- 
tween the substrate and the gate and is separated from 
them by the interface barriers with the resistances R\ 
and R 2l respectively. We consider the stationary situa- 
tion and assume that the Thouless energy corresponding 
to diffusion across the film E^ = D/d 2 f , well exceeds 
all the effective temperatures. Then current across the 
interface separating the substrate and film is 



/=^y'de{^ e (e)- J P|(e)-^(e)+^(e)} ! (71) 



crete example. In the equilibrium F e = tanhe/2T, and rgnt f n ows F e/h (e) = xF e/h (e) + (1 



T e = T cp = T DOS = T Ah = T. Then, we reproduce the 
well-known results for the DOS, MT and AL corrections 
to the conductivity. 

Note that calculation of the DOS, MT and AL cor- 
rections in the equilibrium within the Kcldysh formalism 
was done in Ref. 0. In their derivation of DOS correction 
they have missed the last four terms from Eq. (pq), that 
give important contribution in the final result. However, 
after some canceling mistakes they surprisingly obtained 
the correct final result. 



and a similar equation holds for the interface between 
the film and the gate [here the subscript S denotes the 
substrate]. From the continuity equation for the cur- 

x)Fa (e), where 



F, 



e/h 



R 2 /(Ri + R2). Here F e s /h {e) = tanh (e/(2T x )) and 

tanh ((e =F eV G )/(2T 2 )) denote distributions it the 
substrate and it the gate. Then the gauge invariant dis- 
tribution in the film assumes the form 



F e / h (e) = xta,nh 



e ± (1 - x)eV G 



+(1 — x) tanh 



2T X 

e =F xeVc 



2T 2 



(72) 
(73) 
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FIG. 2: (Color online) The upper figure shows two regions in 
which t G l has different signs, for fixed x — 0.3. In the blue 
region it is greater than zero, and in the white region it is 
smaller than zero. The lower figure shows Tq L dependence on 
Ti, for fixed T 2 and x = 0.3. 



in the case of very resistive interfaces, i.e. when the re- 
sistance of the film can be neglected: R to t ~ R\ + R 2 - 
Next, we calculate parameters appearing in the GL-like 
action, Eq. (|3l|), and demonstrate that they are strongly 
drive dependent. 



A. Ginzburg-Landau relaxation time 



In this subsection we analyze GL relaxation rate un- 
der FFE conditions (given by Eq. (pq)) for the set-up 
shown in Fig. EL The GL relaxation rate is very im- 
portant parameter, since many quantities influenced by 
superconducting fluctuations are singular functions of it. 
One example are fluctuation induced corrections to the 
conductivity, that are analyzed in the previous section. 
While in the equilibrium {T\ = T 2 and V G = Q) and close 
to the transition, it behaves as (|42|), far from equilibrium 
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FIG. 3: (Color online) GL relaxation rate dependence on the 
gate voltage for T\ — T2 = T and x = 0.5. One distinguishes 
quadratic and linear dependence on the gate voltage around 
the transition for T = T c and T — 0.95T C , respectively. 



we find 
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where ^f(z) is the digamma function defined as ^(z) = 
r'(z)/r(z), where T(z) is the gamma function. 7 is the 
Euler con stant and definitions for z cp and F R are given 
in Secjlll[ Note that the theory presented in the previous 
chapters is valid only above (r GL > 0), and very close to 
the transition. However, the obtained expression for GL 
relaxation rate might be valid also below the transition. 
Also, all the expressions are valid for sufficiently small 
voltage drop. The system can be driven from equilib- 
rium due to finite difference of the gate and the substrate 
temperatures and/or due to gate voltage. First, we start 
with the zero voltage case. Then the general expressions 
( Rfm ) simplify significantly, and we obtain fl^ 1 = 0, 
z = 1 and 



-1 8 / T x T c 

r GL = ~ T e X log — - log — 
7T V J 2 J2 



(77) 



Simple analysis shows that t^ 1 is negative (positive) 
when both temperatures T± and T 2 are smaller (greater) 
than the critical temperature T c and can be either posi- 
tive or negative when one of the temperatures is greater 
and another is smaller than T c , see Fig. pi Looking at 
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of parameters, the situation can be more complicated, as 
it is shown in Fig. [|. 



B. Corrections to the electrical conductivity 

In this subsection we examine fluctuation induced cor- 
rections to the in-plane conductivity for the setup in 
Fig. |[ We start with DOS correction, Eq. (p8|). We 
find that characteristic energy scale T oos , Eq. (p(J), pa- 
rameterizing DOS correction reads: 



FIG. 4: (Color online) GL relaxation rate dependence on 
the gate voltage for x — 0.3, T2 = 2T C and T\ — 



T. 



((x-l)/x) T (l/x) 



0.2T C . 



lower part of Fig. g, one notices that GL relaxation rate 
can take rather different values than in the equilibrium. 
The dotted line denotes the equilibrium situation. In the 
equilibrium, we reproduce the (E2) when the system is 
close to the transition (T — T c -C T c ). Next, we switch 
on the gate voltage. Since the expression for GL relax- 
ation rate is rather complicated, we first analyze the case 
Ti=T 2 = T. 



(78) 



where £(2;) is Riemann zeta function. Then, one sees that 
for T = T c and V G = the system is at the transition. 
It moves away from the transition by increasing the gate 
voltage, as it is shown in Fig.pl The GL relaxation rate 
increases quadratically with V G and assumes the form 
t~2 = UV 2 x(l - a;)C(3)/(T c 7r 3 ). On the other hand, if 
for zero gate voltage the system is below, but close to 
the transition {\t^(V = 0,T)| < T < T c ), then at 
some critical finite voltage V c it will be at the transition: 



Then, for V G ; 
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In this case, one obtains linear behavior in V — V c 
around the transition, for the fixed temperature: Tq^ ~ 

/log (^) x(l -x) [^C(3) - logf }{V-V C ). This sit- 
uation is illustrated in Fig. ||. However, for some choice 
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Now, we have the analytic expression for 6a uos as a func- 
tion of T\, T 2 , V G and x. In Fig. ||a we plot dependence 
on the gate and substrate temperatures for zero voltage 
and x = 0.5. In Fig. Sb we plot dependence on gate 
voltage for T\ = T 2 and x = 0.5. 

Next, we analyze AL correction Eq. (|70|). We find that 
the characteristic temperatures T a and T},, Eqs. ( |S8| , |S9|) , 
for the given setup become: 
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(82) 



Then, we have the analytic form of the AL correction. 
Its dependence on system parameters is illustrated in 
Figs, pb, pki. All the temperatures that appear in 5<r MT 
are already calculated. We plot just the first term in 
Eq. (p3|) , the so-called anomalous part of the MT correc- 
tion, in Fig. ||c and Fig. |y, since the second one is equal to 
Sa-Dos. The dephasing rate is taken to be tV 1 = 10~ 3 T C . 
Comparing all the corrections to the conductivity, one 
sees that the most important one close to the transition 
is the MT correction. 



r 



VI. CONCLUSIONS AND DISCUSSION 

We have derived GL-like theory valid close to the tran- 
sition into the superconducting state under FFE condi- 



tions. We considered stationary situation and electro- 
magnetic fields slowly varying (with respect to the su- 
perconducting order parameter) in space. We found that 
the parameters appearing in the GL type action are func- 
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FIG. 5: (Color online) DOS, AL and MT corrections to the in-plane electric conductivity as a function of the gate and the 
substrate temperatures for zero gate voltage, x — 0.5 and tT 1 — 10 T c are shown in Figs, a, c, e, respectively. DOS, AL and 
MT corrections to electric conductivity as a function of gate voltage for Ti = T2 = T, x — 0.5 and tT 1 — 10 _3 T C are shown in 
Figs, b, d, f, respectively. 



tionals of electron and hole noncquilibrium distribution 
functions. Close to an equilibrium, we reproduced known 
results and showed that the time-depenclept| GL equa- 
tion, that is frequently found in liter aturecHlZI is correct, 
contrary to the findings in Reffl. 

Utilizing the theory, we studied the influence of fluc- 
tuations on the electrical conductivity, in FFE situa- 
tion. We demonstrated that different fluctuation induced 
corrections are controlled by different effective tempera- 
tures. These temperatures are drive dependent and carry 
information about temperatures and electrochemical po- 
tentials of the reservoirs that are in contact with the sys- 



tem. We proposed the experimental setup where our pre- 
dictions could be tested, see Fig. [j]. 

The approach developed in the present paper allows 
to analytically treat many other important questions, 
for example, the influence of superconducting fluctua- 
tions on the thermal conductivity under FFE conditions. 
This question, close to an equilibrium, was a contsover- 
sial and puzzling issue for a long time (see Ref.li3 and 
references therein). Hopefully, the final solution is that 
singular contributions of the the DOS and MT correc- 
tion to the thermal conductiidty cancel each other, while 
the AL contribution is finite.Ea Therefore, experimentally 



11 



detected structure in the thermal conductivityE3c3, that 
previously was believed to be explained by these correc- 
tions, needs a new explanation. However, since we have 
found that different corrections to the electrical conduc- 
tivity are characterized by different effective tempera- 



tures under FFE, it is likely that this is the case also 
with thermal conductivity. Then, the MT and DOS cor- 
rections do not cancel each other but instead produce 
nontrivial result, that might explain the experiments. 
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Appendix A: Ginzburg-Landau action 



In this appendix we present a detailed derivation of the nonequilibrium GL-like action Eq. ([3l]) valid in the normal 
state, but very close to the superconducting transition. 

We start from Eq. (fM)]). In the following, we consider stationary distribution functions and stationary electromag- 
netic fields, slowly varying in space with respect to the fluctuating superconducting order parameter A . Then, 



(w £1 , ea (qiX 3)6t (qa)>»s„2 
(w 6l ,e 2 (qi)w* 3>e4 (q2))iS„ 2 
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2(2^) d + 2 8{e 1 - e 4 )S(e 3 - e 2 ) < 5(q 1 + q 2 ) 
irv ' D(qi - 2eA^) 2 + i(e 2 + e 4 ) ' 



(Al) 
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where d denotes the dimension. The average values of all other two field combinations give zero contribution. Then, 
after some algebra we obtain 



S'gl = 2v 



du! dq 



A£_Z^ A£ + + A£_L-/Al + + A£_LJ?AI + 



(A3) 



2ir (2ir) d 
where Aac ]=f = Aj;(^q, T w )- The general formula for the retarded part of superconductive fluctuation propagator is 
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where e_ — e - u/2, Ff = (Fj, 1 - Ff.J/2, and (L A 1 (q, u))* = X iJ 1 (q,cj). In the following we show that L^ can be 
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±ed 



written in the form of Eq. (pj) . Taking into account that gauge invariant distribution function is Ft = F^_ ^ and 
by replacing s = e_ + e<\>K.- onc obtains 

1 fWu pR _ pR « D pR 
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2e 



fU) D 

de 



pR_pR 



D(q - 2eA^) 2 - 2ie - iuj + 2iecj>' 
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Here we have taken into account that Debye frequency ujd 3> u/2 — e<ff^. Also we have added and subtracted the 
term J"° de (P R — -FqM /(2e) from Eq. (A4). Let us concentrate on last two terms in Eq. ( JAg ) and denote their 



sum as l R l : 



i = Z» [°° d J (F e R -F R )[D(q-2eA& 
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The value of the first term in Eq. (A.6) is determined by the poles of the function F R with positive imaginary part. 
Then, very close to the superconductor-metal transition it becomes 



# . f p(, - *as« - 1, + 2^a £ d,(Ml - Of., 



(A7) 



12 



since ui/2 + e0^, as well as D(q — 2eA^) 2 , is much smaller than any relevant scale of the distribution functions. 
Taking into account that f" D detanh (e/2T)/e « 2\og(4:UJr)e 1 /2irT) for oj£> ;2> T where 7 is Euler constant, T c = 

2w£>e 7 ~ 1 / A / 7r an( l Sokhotsky's formula (e + iO)^ 1 = —inS(e) + V{e^ 1 ) one obtains: 



Ln l = 



tanh ( 2^- 



2e 



[--D(q - 2eA£) 2 + iw - 2ie^] + 



TTdF R 
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pi? pi?. 



'4 1 - (A8) 



After introducing T e , $7 and t gl as given by Eqs.( [74[]37| ), on e arrives at Eq. (|34j). 



Note that the term ~ A d *A ci van i snes . This is expected 
property of the action (see the explanation in the main 
text) . Taking into account that dt in Eq.Q is just a symbol 
standing instead of a matrix in the discrete time space, 
and that sum of the retarded and the advanced Green's 
functions taken at the same time vanishes, we find that 
term ~ A c/ *A d is zero, while 



L K ± = 1 I de 



1 _ A I ph F e 



rph, pe 



D(q-2eA^) 2 - 2ie 
Very close to the transition, it reduces to 
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Note that after applying the Sokhotsky's formula, the 



term/ de(F^F^ +F^ e Ff e _)/e gives zero contribution, 
since the function under the integral is the odd function 

of e. 



In the above calculation we used that close to the 
transition cu/2 — ecjfj^, as well as Z3(q — 2eA^) 2 , is 
much smaller than any relevant scale of the distribu- 
tion functions. Let us demonstrate the importance of 
the mentioned combination of the momenta/frequancy 
and gauge fields by considering simple example. Let 
us consider an equilibrium situation and denote the or- 
der parameter by Ao(w,q) in this case. It is satisfied 



Dq 2 



scalar potential 

A = Aq cxp (—2ie<fit) 



(tgl) ' < T. Next we turn on a constant 
Then, the order parameter becomes 



i.e. A(w,q) = A (w — 2ecj), q). 
Then, uj — 2ecf> ~ (tg2) _1 ^C T. Similarly, it can be shown 
that the combination D(q — 2eA^) 2 has to be compared 
with characteristic energies of the distribution function. 
That is why it is necessary to introduce the gauge invari- 
ant distribution functions during the calculation, as wc 
did above. 
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